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The arctic circle theorem of Jockusch, Propp, and Shor asserts that 
uniformly random domino tilings of an Aztec diamond of high order 
^ are frozen with asymptotically high probability outside the "arctic cir- 

cle" inscribed within the diamond. A similar arctic circle phenomenon 
y—i has been observed in the limiting behavior of random square Young 

tableaux. In this paper, we show that random domino tilings of the 
Aztec diamond are asymptotically related to random square Young 
tableaux in a more refined sense that looks also at the behavior in- 
^ J side the arctic circle. This is done by giving a new derivation of the 

limiting shape of the height function of a random domino tiling of the 
Aztec diamond that uses the large-deviation techniques developed for 
O the square Young tableaux problem in a previous paper by Pittel and 

^ the author. The solution of the variational problem that arises for 

domino tilings is almost identical to the solution for the case of square 
Young tableaux by Pittel and the author. The analytic techniques 
used to solve the variational problem provide a systematic, guess-free 
approach for solving problems of this type which have appeared in a 
number of related combinatorial probability models. 
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matrix, Aztec diamond, arctic circle, large deviations, variational problem, 
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1 Introduction 

1.1 Domino tilings and the arctic circle theorem 

A domino in is a Z^-translate of either of the two sets [0, 1] x [0, 2] or 
[0, 2] X [0, 1]. If 5* C is a region comprised of a union of Z^-translates of 
[0, 1]^, a domino tiling of 5* is a representation of 5 as a union of dominoes 
with pairwise disjoint interiors. Domino tihngs, or equivalently the dimer 
model on a square lattice, are an extensively studied and well-understood 
lattice model in statistical physics and combinatorics. Their rigorous analysis 
dates back to Kasteleyn [18] and Temperley-Fisher [26], who independently 
derived the formula 



for the number of domino tilings of an n x m rectangular region. About 
thirty years later, a different family of regions was found to have a much 
simpler formula for the number of its domino tilings: if we define the Aztec 
diamond of order n to be the set 



domino tilings. This can be proved by induction in several ways, but is 
perhaps best understood via a connection to alternating sign matrices. 

One of the best-known results on domino tilings is the arctic circle the- 
orem due to Jockusch, Propp and Shor [H] , which describes the asymptotic 
behavior of uniformly random domino tilings of the Aztec diamond. Roughly, 
the theorem states that the so-called polar regions, which are the four con- 
tiguous regions adjacent to the four corners of the Aztec diamond in which 
the tiling behaves in a predictable brickwork pattern, cover a region that is 
approximately equal to the area that lies outside the circle inscribed in the 
diamond. See Figure |2| where the outline of the so-called "arctic" circle can 
be clearly discerned. The precise statement is the following. 






Elkies, Kuperberg, Larsen and Propp [U] proved that 



AD„ has exactly 
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(a) (b) 
Figure 1: The Aztec diamond of order 3 and one of its 64 tilings by dominoes. 



Theorem 1 (The arctic circle theorem Fix e > 0. For each n, consider 
a uniformly random domino tiling of ADn scaled by a factor 1/n in each axis 
to fit into the limiting diamond 

AD^:={\x\ + \y\<l}, 

and let P° C n~^ADn he the image of the polar regions of the random tiling 
under this scaling transformation. Then as n ^ oo the event that 

G AD^:x'' + y^ > ^ + n {n~' AD„) 

C P: C !^ix,y) eAD^:x^ + y^>^-e^ 

holds with probability that tends to 1. 

In later work, Cohn, Elkies and Propp [2] derived more detailed asymp- 
totic information about the behavior of random domino tilings of the Aztec 
diamond, that gives a quantitative description of the behavior of the tiling 
inside the arctic circle. They proved two main results (which are roughly 
equivalent, if some technicalities are ignored), concerning the placement 
probabilities (the probabilities to observe a given type of domino in a given 
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Figure 2: The arctic circle theorem: in a random domino tihng of AD50, the 
circle-hke shape is clearly visible. Here, dominoes are colored according to 
their type and parity. 



position in the diamond) and the height function of the tiling (which, 
roughly speaking, encodes a weighted counting of the number of dominoes 
of different types encountered while travelling from a fixed place to a given 
position in the diamond - see Section [6] for the precise definition). 

A main goal of this paper is to give a new proof of the Cohn-Elkies- 
Propp limit shape theorem for the height function of a uniformly random 
domino tiling of the Aztec diamond - see Theorem 12 in Section [6j Our 
proof is based on a large deviations analysis, and so gives some information 
that the proof in ^ (which is based on generating functions) does not: a 
large deviation principle for the height function. Perhaps more importantly, 
it highlights a surprising connection between the domino tilings model and 
another, seemingly unrelated, combinatorial probability model, namely that 
of random square Young tableaux. 
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Figure 3: A square Young tableau of order 5 (shown in the "French" coordi- 
nate system), and the wall whose construction the tableau encodes at various 
stages of its construction. 



1.2 Random square Young tableaux 

Recall that a square (standard) Young tableau of order n is an array (tij)i^j=i 
of integers whose entries consist of the integers 1,2,..., n^, each one appear- 
ing exactly once, and such that each row and column are arranged in increas- 
ing order. One can think of a square Young tableau as encoding a sequence 
of instructions for constructing an nx n wall of square bricks leaning against 
the X— and y— axes by laying bricks sequentially, where the rule is that each 
brick can be placed only in a position which is supported from below and 
from the left by existing bricks or by the axes. In this interpretation, the 
number tij represents the time at which a brick was added in position 
see Figure |3} The number of square Young tableaux of order n is known (via 
the hook-length formula of Frame-Thrall-Robinson) to be 



n:,=i(^+j-i)' 

In [22], Boris Pittel and the author solved the problem of finding the 
limiting growth profile, or limit shape, of a randomly chosen square 
Young tableau of high order. In other words, the question is to find the 
growth profile of the square wall "constructed in the most random way". 
This can be expressed either in terms of the limit in probability L{x,y) of 
the scaled tableau entries n~Hij, where {x,y) G [0,1]^ and i = i{n) and 
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Figure 4: The limiting growth profile of a random square Young tableau and 
the profile of a randomly sampled tableau of order 100. The curves shown 
correspond to (scaled) times t = i/lO, z = 1, 2, . . . , 9. 



j = j{n) are some sequences such that i/n ^ x and j /n y as n —>■ oo; or 
alternatively in terms of the limiting shape of the family of scaled "sublevel 
sets" {n~^{i,j) : tij < a-in?} for each a G (0, 1) (which in the "wall-building" 
metaphor represents the shape of the wall at various times, and thus can be 
thought of as encoding the growth profile of the wall). Figure |4] shows the 
limiting growth profile found by Pittel and Romik and the corresponding 
profile of a randomly sampled square Young tableau of order 100. 

For the precise definition of the limiting growth profile, see [22]. Here, we 
mention only the following fact which will be needed in the next subsection: 
If L : [0, 1] X [0, 1] [0, 1] is the limit shape function mentioned above, then 
its values along the boundary of the square are given by 

L(0,t) = L{t,0) = ^'"(y^ , (0<t<l), (1) 



L(l,t) = L(t,l) = ^ + *\ (0<t<l). (2) 

Also note that according to the limit shape theorem, the convergence of 
to L{i/n, j /n) as n oo is uniform in i and j (this follows easily 
from monotonicity considerations). 
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1.3 An arctic circle theorem for square Young tableaux 

While it is not immediately apparent from the description of this limit 
shape result, it follows as a simple corollary of it that random square Young 
tableaux also exhibit an "arctic circle" -type phenomenon. That is, there is 
an equivalent way of visualizing the random tableau in which a spatial phase 
transition can be seen occurring along a circular boundary, where outside 
the circle the behavior is asymptotically deterministic (the "frozen" phase) 
and inside the circle the behavior is essentially random (the "disordered" or 
"temperate" phase). This fact, overlooked at the time of publication of the 
paper [22], was observed shortly afterwards by Benedek Valko [27j. In fact, 
deducing the arctic circle result is easy and requires only the facts Q, (g 
mentioned above, which contain only a small part of the information of the 
limit shape. 

To see how the arctic circle appears, we consider a different encoding of 
the information contained in the tableau via a system of particles on the 
integer lattice Z. In this encoding we have n particles numbered 1, 2, . . . , n, 
where initially, each particle with index k is in position k. The particles are 
constrained to remain in the interval [1, 2n\. At discrete time steps, particles 
jump one step to the right, provided that the space to their right is empty 
(and provided that they do not leave the interval [1, 2n]). At each time step, 
exactly one particle jumps. 

It is easy to see that after exactly rt^ steps, the system will terminate 
when it reaches the state in which each particle k is in position n + k, and 
no further jumps can take place. We call the instructions for evolving the 
system of particles from start to finish a jump sequence. We can now add a 
probabilistic element to this combinatorial model by considering the uniform 
probability measure on the set of all jump sequences of order n, and name the 
resulting probability model the jump process of order n. But in fact, this 
is nothing more than a thinly disguised version of the random square Young 
tableaux model, since jump sequences are in a simple bijection with square 
Young tableaux: given a square tableau, think of the sequence of numbers 
in row k of the tableau as representing the sequence of times during which 
particle n + 1 — k jumps to the right. This is illustrated in Fig. |5| We leave 
to the reader the easy verification that this gives the desired bijection. 

With these definitions, it is now natural to consider the asymptotic be- 
havior of this system of particles as tt, — cxd. Figure |6] shows the result for 
a simulated system with n = 40. Here we see a circle-like shape appearing 
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Figure 5: The bijection between square Young tableaux and jump sequences: 
Each row in the tableau encodes the sequence of times at which a given 
particle jumps. As an example, the highlighted trajectory on the right-hand 
side corresponds to the highlighted row on the left-hand side. 



again. To formulate precisely what is happening, given a jump process of 
order n, for each I < k < 2n, let r~(/c) and t+(A;) denote respectively the 
first and last times at which a particle k jumped from or to position k. Define 
the frozen time-period in position k to be the union of the two intervals 

[0,T-ik)]U[r^{k),n']. 
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Figure 6: A jump process with 40 particles. 



Theorem 2 (The arctic circle theorem for random square Young tableaux) . 
Fix any e > 0. Denote 



holds with probability that tends to 1. In other words, if the space-time di- 
agram of the trajectories in a random jump process is mapped to the unit 
square [0, 1] x [0, 1] by scaling the time axis by a factor l/n^ and scaling the 
position axis by a factor of l/2n, then for large n the frozen time-periods 
will occupy approximately the part of the space-time diagram that lies in the 
complement of the disc 



V±{x) = - ± Va;(l -x). 



As n 



oo, the event 





inscribed in the square. 
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Proof. First, note the following simple observations that express the times 
T~{k) and r+(/c) in terms of the Young tableau {ti,j)^^j=i- 

(i) For 1 < A; < n we have T~{k) = tn+i-k,i- 

(ii) For n + 1 < A; < 2ri we have T~{k) = ti^k-n- 

(iii) For 1 < A; < n we have t+(/c) = tn,k- 

(iv) For n + 1 < k < 2n we have T^{k) = t2n+i-k,n- 

For example, the first statement is based on the fact that when 1 < k < n, 
the time T~{k) is simply the first time at which the particle starting at 
position k (which corresponds to row n + 1 — in the tableau) jumps. The 
three remaining cases are equally simple and may be easily verified by the 
reader. 

Combining these observations with ([T]) and (|2| and the limit shape the- 
orem, we now see that after scaling the times T~{k) and T^{k) by a factor 
of n~^, we get quantities that converge in the limit, uniformly in k, to val- 
ues determined by the appropriate substitution of boundary values in the 
limit shape function L{x., y). For example, to deal with case (i) above, when 
1 < k < n, using Q we have that 

/ k — 1' 

^~'^'^n{k) = n~Hn+l-k,l ^ L 0, 1 - 



n 



1 - 



2 2 

uniformly in 1 < k < n. Similarly, the other three cases each imply that 
n~'^T^{k) is uniformly close to ip±{k/2n) in the appropriate range of values 
of k; we omit the details. Combining these four cases gives exactly that the 
event in Theorem [2] holds with asymptotically high probability as n — oo. 

□ 



1.4 Similarity of the models and the analytic technique 

Apart from giving a new proof of the limit shape theorem of Cohn, Elkies 
and Propp, another main goal of this paper is to show that the two mod- 
els described in the preceding sections (random domino tilings of the Aztec 
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diamond and random square Young tableaux) exhibit similar behavior on 
a more detailed level than that of the mere appearance of the arctic circle, 
and that in fact they are almost equivalent in an asymptotic sense. Our new 
proof of the limit shape theorem for the height function will use the same 
techniques developed in [22] for the case of random square Young tableaux: 
we first derive a large deviations principle, not for domino tilings but for a 
related model of random alternating sign matrices, then solve the re- 
sulting problem in the calculus of variations using an analysis that parallels, 
to a remarkable (and, in our opinion, rather surprising) level of similarity, 
the analysis of the variational problem in [22] • The resulting formulas for the 
solution of the variational problem are almost identical to the formulas for 
the limiting growth profile of random square Young tableaux. Up to some 
trivial scaling factors related to the choice of coordinate system, the formulas 
for the two limit shapes can be written in such a way that the only difference 
between them is a single minus sign. 

Another important aspect of our results lies not in the results themselves 
but in the techniques used. We use the methods first presented in [22] to 
solve another variational problem belonging to a class of problems previously 
thought to be difficult to analyze, due to a lack of a systematic framework 
that enables one to derive the solution in a relatively mechanical way (as 
opposed to having to guess it using some deep analytic insight) and then 
rigorously verify its claimed extremal properties. This justifies to some extent 
the claim from [22] that the analytic techniques of that paper provide a 
systematic approach for dealing with such problems, which seem to appear 
frequently in the analysis of combinatorial probability models (see [3l [191 1221 
I2H1 [23 ) , and are also strongly related to classical variational problems arising 
in electrostatics and in random matrix theory. 

The rest of the paper is organized as follows. In Section [2] we recall some 
facts about alternating sign matrices, and study the problem of finding the 
limiting height matrix of an alternating sign matrix chosen randomly ac- 
cording to domino measure, which is a natural (non-uniform) probability 
measure on the set of alternating sign matrices of order n. In Section [3] we 
derive a large deviation principle for this model. This problem is solved in 
Section |4} In Section [5] we prove a limiting shape theorem for the height ma- 
trix of an alternating sign matrix chosen according to domino measure. In 
Section[6]we deduce from the previous results the Cohn-Elkies-Propp limiting 
shape theorem for the height function of uniformly random domino tilings of 
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Figure 7: (a) An ASM of order 6; (b) its height matrix. 



the Aztec diamond. Section [7] has some final remarks, including a discussion 
on the potential applicability of our methods to attack the well-known open 
problem of the limit shape of uniformly random alternating sign matrices. 

2 Alternating sign matrices 

An alternating sign matrix (often abbreviated as ASM) of order n is 
an n X n matrix with entries in {0, —1, 1} such that in every row and every 
column, the sum of the entries is 1 and the non-zero numbers appear with al- 
ternating signs. See Fig. [7]^a) for an example. Alternating sign matrices were 
first defined and studied in the early 1980's by David Robbins and Howard 
Rumsey in connection with their study [25j of Charles Dodgson's conden- 
sation method for computing determinants and of the X- determinant, a 
natural generalization of the determinant that arises from the condensation 
algorithm. Later, Robbins, Rumsey and William Mills published several in- 
triguing theorems and conjectures about them [SD], tying them to the study 
of plane partitions and leading to many later interesting developments, some 
of which are described, e.g., in [HEl]. 

Denote by An the set of ASM's of order n. For a matrix M e An, denote 
by N^{M) the number of its entries equal to 1. An important formula proved 
by Mills, Robbins and Rumsey states that 

J2 2^+W = 2("^'). (3) 

M e^„ 

This is sometimes referred to as the "2-enumeration" of ASM's. The reader 
may note that the right-hand side is equal to the number of domino tilings of 
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AD„ mentioned at the beginning of the introduction; indeed, a combinatorial 
explanation for (|3| in terms of domino tilings was found by Elkies, Kuperberg, 
Larsen and Propp |j6j. In Section |6] we will say more about this connection 
and how to make use of it, but for now, we rewrite ^ more probabilistically 

as 

2-("J') J2 ^^^^'''^ = 1' 

and consider this as the basis for defining a probability measure on An, 
which we call domino measure (thus named since it is closely related to 
the uniform measure on domino tilings of AD„; see Section [6]), given by the 
expression 

PSom(M) = 2^+(^^^)~("J'), {Me An). 

Our first goal will be to study the asymptotic behavior of large random 
ASM's chosen according to domino measure, and specifically the limit shape 
of their height matrix. The height matrix of an ASM M = {iTLij)^j^i G An 
is defined to be the new matrix H{M) = {hij)'^j^Q of order (n + 1) x (n + 1) 
whose entries are given by 



p<i q<j 



Tflpq. 



The matrix H{M) is also sometimes referred to as the corner sum matrix 
of M. It satisfies the following conditions: 

^o,fc = hkfi = for all < k < n, (HI) 

hn,k = hk^n = k for all < A; < n, (H2) 
< /ij+ij — /ijj, — hj^i < 1 for all < z < < j < ra. (H3) 

See Fig. [7](b) for an example. (In fact, it is not too difficult to see that the 
correspondence M — > H{M) defines a bijection between the set of ASM's 



of order n and the set of matrices satisfying conditions (II1)-(II3) - see 



[25l Lemma 1] - but we will not need this fact here). In particular, the 



"Lipschitz"-type condition (H3) means that the height matrix can be thought 
of as a discrete version of a two-dimensional surface, and is therefore a natural 
candidate for which to try and prove a limit shape result. 

The basis for our analysis of Pp^^j-random ASM's is a formula which will 
give the probability distribution (under the measure Pdom) of ^^e k-th row 
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of the height matrix, for each 1 < k < n. To describe this, first, as usual, 
denote the Vandermonde function by 

A{ui, . . . ,Um) = Y\. i^j-Ui). 

l<i<j<m 

Second, for an ASM M e An and some 1 < A; < n, let Xfc(l) < Xk{2) < 
. . . < Xjs{k) be the unique ascents of the A;-th row of the height matrix 
H{M), namely those column indices such that 

hk,Xk{i) " hk,Xkii)-i = 1) {i = 1,2, . . . ,k). 



Note that the conditions (H1)-(H3) guarantee that the ascents exist, that 
there are exactly k of them, and that the original k-th row of H[M) can be 
recovered from them. 

Theorem 3. // integers 1 < Xi < X2 < • . • < < n are given, and if 
Vi < 1/2 < ■ ■ ■ < Vn-k o^re the numbers in {1, 2, . . . , n} \ {xi, . . . , Xk] arranged 
in increasing order, then, in the notation above, we have 



DOM 



M eAn-. {Xk{l), Xk{k)) = {xi, ...,Xk) 

{k + l\ (n-k + l\ . , ^ . ^ 

2l 2 j2l 2 ) A{xi,...,Xk)A{yi,...,yn-k) 



2(T) A{l,2,...,k)A{l,2,...,n-k) 



■ (4) 



To prove Theorem [3} we use another well-known combinatorial bijection 
relating ASM's to monotone triangles . A monotone triangle of order n is 
a triangular array {tij)i<i<n,i<j<i of integers satisfying the inequalities 

tij < tij+i, tij < ti^ij < tij+i {2 < i < n,l < j < i - 1). 

A complete monotone triangle of order n is a monotone triangle whose 
bottom row consists of the numbers (1,2, . . . ,n). It is well-known that al- 
ternating sign matrices of order n are in bijection with complete monotone 
triangles of order n. In our terminology, the bijection assigns to an ASM 
M = {mij)^j^i the monotone triangle 

^ — i^i,j)^<i<n,l<j<i = V^ASM^CMtIA^ ) 

whose fc-th row {tk,j)i<j<k consists for each 1 < k < n of the ascents of 
the k-th row of the height matrix H{M), arranged in increasing order. See 
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Figure 8: (a) The complete monotone triangle corresponding to the ASM in 
Figure [7| (b) its dual, shown "standing on its head" . 



Figure jSj^a) for an example. More explicitly, it is easy to check that this 
means that an index j will be present in the k-th row of T if and only if 

k 

J2 m^J = 1 

i=l 

holds. 

Another notion that will prove useful is that of the dual of a complete 
monotone triangle. If T is a complete monotone triangle of order n, and M 
is the ASM in An such that T = v'asm^cmt(-^), then the dual T* of T is 
the complete monotone triangle of order n that corresponds via the same 
bijection to the matrix W, defined as the vertical reflection of M, i.e., the 
matrix such that Wij = rrin+i-ij for all i,j (clearly it, too, is an ASM). See 
Figure |8](b), where the dual triangle is drawn reflected vertically. 

The following simple observation describes more explicitly the connection 
between a monotone triangle and its dual. 

Lemma 4. IfT= {tij)i<i<n,i<j<i is a complete monotone triangle of order 
n, then for each 1 < k < n — 1, the {n — k)-th row of the dual triangle T* 
consists of the numbers in the complement 

{1, 2, . . . , n} \ {tfc,i, tfc,2, • • • , tk,k} 
of the k-th row of T , arranged in increasing order. 

Proof. Let M = {mij)ij G An be such that T = v^asm^cmt(^)- As 
mentioned above, 1 < j < n appears in the fc-th row of T if and only if 
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^i=i'^i.j = 1- Similarly, from the definition of T* we see that j appears 
in the {n — k)-th row of T* if and only if X^ILfc+i^^iJ ~ ^- from the 

definition of an alternating sign matrix, one and only one of these conditions 
must hold. □ 

As the last step in the preparation for proving Theorem |3} we note that 
if M G An and T = (y9ASM^CMT(Af ), then it is easy to see that A^+(M), the 
number of +1 entries in M, can be expressed in terms of T as the num- 
ber of entries tij in T that do not appear in the preceding row (including, 
vacuously, the singleton element in the top row). We denote this quantity 
also by N^{T); note that it is defined more generally also for non-complete 
monotone triangles. We furthermore recall the following formula proved by 
Mills, Robbins and Rumsey in [20l Th. 2] (see also P Eq. (7), Section 4], 
and see [11] for a recent alternative proof and some generalizations): 

Lemma 5. If k > 1 and xi < X2 < ■ ■ ■ < Xk are integers, then the sum of 
2^+(^) oyer all monotone triangles T of order k with bottom row (xi, . . . , x^) 
is equal to 

2(T) 17 ^l^, 

l<i<j<k 

Proof of Theorem^ Denote by Tn{xi, . . . , Xk) the set of complete monotone 
triangles of order n whose /c-th row is equal to (xi, . . . , x^). From the remarks 
above, it follows that the left-hand side of ^ is equal to 



Ter„(a;i,...,a;fc) 

In addition, for a monotone triangle T G 7^(xi, . . . , x^), define Ttop and 
^bottom as the two monotone triangles, of orders k and n — k, respectively, 
where T^p is comprised of the top k rows of T, and Tbottom is comprised of 
the top n ~ k rows of the dual triangle T* . From Lemma |4[ it follows that 
the correspondence 

T (Ttop; ^bottom) 

defines a bijection between 7^(xi, . . . ,Xk) and the cartesian product Ax B, 
where A is the set of monotone triangles with bottom row (xi, . . . , Xk) and 
B is the set of monotone triangles with bottom row (j/i, . . . ,yn-k) (in the 
notation of Theorem [3]). This corresponence furthermore has the property 
that 

iV+(r) = iV+(Ttop) + iV+(Tbottom) 
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(since N^{Ttop) counts the number of +1 entries in tlie first k rows of tlie ASM 
corresponding to T, wfiereas A^+(Tbottom) counts tfie number +l's in tlie last 
n — k rows), or equivalently that 2^+^-^^ = 2^+(-^*°p)2^+(-^''°"°™^ Combining 
these last observations, we get that the left-hand side of ^ is equal to 

2~i"t^) ^ ^ 2-'^+(-^'°p) ^ ^ 2''^+'- -^bottom) 

which by Lemma [s] is equal exactly to the right-hand side of (|4]) . □ 

We remark that an equivalent version of Theorem |3} phrased in the lan- 
guage of domino tilings and certain so-called zig-zag paths defined in terms 
of them, is proved by Johansson in [T3] (see Proposition 5.14 in that paper 
and eq. (5.16) following it). See also the subsequent papers [Ml [IZ] where 
Johansson proves many interesting results about random domino tilings of 
the Aztec diamond by combining a variant of ^ with ideas from the theory 
of orthogonal polynomials and the theory of determinantal point processes. 



3 A large deviation principle 

We now turn from combinatorics to analysis, with the goal in mind being 
to use Theorem [3] as the starting point for a large deviation analysis of the 
behavior of P^Q^-random ASM's. First, we define the space of functions on 
which our analysis takes place. Fix < y < 1. We wish to understand the 
behavior of the k-th row of the height matrix of a PJ^Qj^-random ASM of 
order n for values of k satisfying k ^ y ■ n, when n is large. 
Define the space of y-admissible functions to be the set 

JFj, = |/ : [0, 1] — s> [0, 1] : / is monotone nondecreasing, 1-Lipschitz, 

and satisfies /(O) = 0, /(I) = y|. 

Define the space of admissible functions as the union of all the ?/-admissible 
function spaces: 

^= U 

ye [0,1] 
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Figure 9: A (6, 3)-admissible sequence u and the corresponding function /„. 



We also define a discrete analogue of the admissible functions. Given 
integers < k < n, a sequence u = (mo,mi, . . . ,Un) of integers is called an 
{n, k) -admissible sequence if it satisfies 

Mo = 0, Un = k, and Wj+i — Ui & {0, 1} for all < i < n — 1, 

Note that (n, /c)-admissible sequences are exactly those that can appear as 
the k-th row of a height matrix H{M) of an ASM M G An- We embed the 
(n, fc)-admissible sequences in the space J^y for y = k/n, in the following way: 
For each (n, fc)- admissible sequence u, define a function : [0, 1] [0, 1] as 
the unique function having the values 

fuij/n)=Uj/n, 0<j<n, 

and on each interval [j/n, {j + l)/n] for < j < n — 1 being defined as 
the linear interpolation of the values on the endpoints of the interval; see 
Figure [9j Clearly, /u is a (fc/?T,)-admissible function. In fact, it is easy to see 
that the admissible functions are precisely the limits of such functions in the 
uniform norm topology. 

With these definitions, we can now formulate the large deviation principle. 
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Theorem 6 (Large deviation principle for P^Qj^j-random ASM's). Let < 

k < n, and let u = (mq, mi, . . . , «„) be an {n, k) -admissible sequence. Let 
H{M)k denote the k-th row of a height matrix H{M). Then 



p: 



DOM 



M eAn-- H{M)k = u 



exp ( - (1 + o(l))n2(/(/„) + eik/n))) , (5) 



where we define 
1 



oiy) 



2 log 2 - 3 



3 
2' 



- r [\og\s-t\f\s){f\t)-l)dsdt, (/e^). 

Jo Jo 



The o(l) error term in ([s]) is uniform over all < k < n and all {n,k)- 
admissible sequences u, as n ^ oo. 

Proof. Let 1 < xi < X2 < ■ ■ ■ < < n he the positions of the k ascents in 
the sequence {uq, Ui, . . . , Un) (in the same sense defined before, namely that 

— = 1); and let 1 < yi < . . . < yn-k < nhe the numbers in the 
complement {1, . . . , n} \ {xi, . . . , Xk] arranged in increasing order. 

By (|4]), we have 



n-' log PSo^ 



Me An-. H{M)k = u 



"-(cr)H"""0-cr))-^ 

l<i<j<k l<i<j<n—k 

+n"^ ^ log(xj - Xi) + n"^ ^ log(?/j - yi) 

l<i<j<k l<i<j<n—k 



(6) 
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We estimate each of the summands. First, we have 

log2--fl--')+o(l). (7) 
n \ n I 



Secondly, the sum n ^ X]i<j<j<A: log(j — i) can be rewritten as 

k-l 

^ \og{] ~i) =n''^^{k- d)\ogd 

'.i<j<k k—l d=l /j\2 k—1 ^ ,\ , _ 

(i=l ^ ^ (i=l ^ ^ 

- ^^^-^\ogk+(^^ j\l~t)\ogtdt + o{l) 



2n2 

1 , 3 /A;^ ^ 



logA;-- - +o(l), (8) 
n J 4 \?T, / 

where the error term o(l) is uniform in as n ^ oo (the estimate for this 
sum is essentially the leading-order asymptotic expansion for the so-called 
Barnes G -function, related also to the hyperfactorial; for more detailed 
asymptotics of these special functions, see ^ Sec. 2.15, p. 135]). Similarly, 
replacing A; by n — we get that 

l<i<j<n-k ^ ^ ^ ^ 

Finally, we estimate the terms in ^ that depend directly on the sequence u. 
The idea is to replace each term n"^ log(a;j — Xi) with an integral of the form 
JJ\og{t — s) f^(t) f^{s) ds dt over a certain region. Observe that for X > 1 
we have the (easily verifiable) identity 

1 rX+l 

\og{v — u) dv du 

^ logX + f hx^ + 1) log ( + X log + 



Jx 



2' ' ^\ °VX-1/ 2 



21 



When X is large, this behaves hke log X + (^) . The integral is also defined 
and finite when X = 1. So we can write 



n 2 loj 

l<i<j<A; 



[^Xj Xi) 



n 



E 



n 



i<«<j<fe 

E 



log(i; — u) dv du + O \ } 

\i<*<i<fc ■' 
logn' 



dvdu 
log(f - u) — — + O 



n 



l<i<j<k "'^i-i 



, /v — u\ dvdu , k(k — l] 
log ( ^— ) — ^ + logn ■ 



n 



+0 



log 77. 



Now observe that j'Jyx) is equal to 1 if (xj — l)/r;, < a; < Xijn for some i, or 
to otherwise; so this last expression can be rewritten as 



U \o^(t-s)Us)Ut)dsdt^ 



2 \n I \ n 



(10) 



where the integral is over the region 

U ' 

l<i<j<™ 



n n 



3 - 1 J 



n n 



The region of integration in (10) can be replaced with the slightly larger 
region 

R= {(s,t) G [0,1] X [0, 1] : s < t}, 
at the cost of an additional error which can be bounded in absolute value by 



1 ry 

dy / \\og{y — x)\dx 

Jy-l/n 



1/n 



log tdt 



O 



logn 



n 



To summarize, after this change in the region of integration and, in addition, 
after symmetrizing the region of integration for convenience, we have shown 
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that 



l<i<j<k 
»1 rl 



JO 



\og\t-s\as)mdsdt+U'^] logn + Or^]. (11) 



Symmetrically, following exactly the same reasoning for the last sum in ^ 
we get the similar estimate 



l<i<j<n—^k 
>1 /.I 



Vi) 



n 



Ol'^'\. (12) 



n 



It remains to plug the estimates ([T]), ([s]), ([9]), (11) and (12) into ([6]), and 
simplify. Denoting y = k/n, and using the integral evaluation 

- I I log It — s\ds dt - 
2 Jo Jo 

this gives that the left-hand side of ^ is equal to 



3 

■4' 



~ i + i^' + i^^ " " ■ " ^) " logy - ^(1 - yf log(l - y) 



[ [ \og\t~s\fi{s)f:,{t)dsdt- [ [ \og\t-s\fi{s)dsdt + o{l) 
Jo Jo Jo Jo 

= -%)-/(/„) + 0(1), (13) 
as claimed. □ 
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4 The variational problem and its solution 

Fix < y < 1. Motivated by Theorem [6} we now turn our attention to the 
problem of minimizing the integral functional /(/) over the appropriate class 
of ?/-admissible functions. In the next section we will show how this implies 
a limit shape result for Pp^^j-random ASM's. 

The precise variational problem that we will solve is the following: 
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Variational Problem 1. For a given < y < 1, find the function f* that 
minimizes I{f) over all functions f E J^y. 

Variational Problem [T] is a variant of a class of variational problems that 
have appeared in several random combinatorial models; see, e.g., [21 113 [221 
l28l [29] . Such problems bear a strong resemblance to classical physical prob- 
lems of finding the distribution of electrostatic charges subject to various con- 
straints in a one-dimensional space, as well as to problems of finding limiting 
eigenvalue distributions in random matrix theory. However, the variational 
problems arising from combinatorial models usually have non-physical con- 
straints that make the analysis trickier. In particular, in several of the works 
cited above, the presence of such constraints required the authors to first 
(rather ingeniously) guess the solution. Once the solution was conjectured, 
it was possible to verify that it is indeed the correct one using fairly standard 
techniques. Cohn, Larsen and Propp, who derived the limit shape of a ran- 
dom boxed plane partition, ask (see Open Question 6.3 in [3]) whether there 
exists a method of solution for their problem that does not require guessing 
the solution. 

In [22], it was argued however that when dealing with such problems, 
it is not necessary to guess the solution, since a well-known formula in the 
theory of singular integral equations for inverting a Hilbert transform on 
a finite interval actually enables mechanically deriving the solution rather 
than guessing it, once certain intuitively plausible assumptions on the form 
of the solution are made. Here, we demonstrate again the use of this more 
systematic approach by using it to solve our variational problem. As an 
added bonus, the solution rather elegantly turns out to be nearly identical to 
the solution of the variational problem for the square Young tableaux case 
(although we see no a priori reasons why this should turn out to be the 
case), and we are able to make use of certain nontrivial computations that 
appeared in [22], which further simplifies the analysis. 

Our goal in the rest of this section will be to prove the following theorem. 
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Theorem 7. Define 
2 



Z{x,y) 



TT 



X — 1/2) arctan 



\-ix-l/2y-{y-l/2y 



1/2-2/ 



+ 1 arctan' 2(x - l/2)(l/2 - ,) 
2 \ /I _ _ 1/2)2 _ _ 1/2)2 

X- 1/2 



(14) 



-(1/2 — y) arctan 



yi - (x - 1/2)2 - (y - 1/2) 
For < ?/ < 1/2, t/ie solution f* to Variational Problem^is given by 







< a; < 



2 ' 2 

2/ 



2 ' 

1+2^^(1-^) 
2 ^ ^ ^ 2 ' 

l+2y/yil-y) 



I + (x, y) ^ < X < 



(15) 



< X < 1. 



X 



For y = 1/2, the solution is given by 
For y > 1/2 the solution is expressed in terms of the solution for 1 — y by 

/; = ^ - /r-r 

Moreover, for all < y < 1 we have 

liFy) = -e{y). 



As a first step, for convenience we reformulate the variational problem 
slightly to bring it to a more symmetric form, by replacing each / G JFj, by 
the function 

g{x) = 2 fix) - X. (16) 

It is easy to check how the class of y-admissible functions and the functional 
/(■) transform under this mapping. The result is the following equivalent 
form of our variational problem. 
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Variational Problem 2. For < y < 1, define the space of functions 
Qy = : [0, 1] [—1, 1] : g{0) = 0, g{l) = 2y — 1, and g is 1-Lipschitz^ 
and the integral functional 



Ji.9) = -j / g{,s)g'{t)\og\s-t\dsdt. 
Jo Jo 

Find the function g* G Qy that minimizes the functional J over all functions 
g e Gy. 



The reader may verify that if / G and g E Qy are related by (16), then 
the integral functionals / and J are related by 

/(/) = - 1 

This implies that the following theorem is an equivalent version of Theorem [7] 

Theorem [T]'. For < y < 1/2, the solution g* to Variational Problem^is 
given by 



-X 



y-x + Z{x,y) 
2y — X 



<x < 

l-2^yy{l-y) 
2 

l+2y/y{l-y) 



l-2y/y{l-y) 



< X < 



l+2y/y{l-y) 



< X < 1, 



where Z{x,y) is defined in (14). For y = 1/2, the solution is given by 
gy^[x) = 0. For y > 1/2 the solution is expressed in terms of the solution 
for 1 — y by g* = —g\_y. Moreover, for all < y < 1 we have 



We now concentrate our efforts on proving Theorem 7' First, in the fol 



lowing lemma we recall some basic facts about the space Qy and the functional 
J. We omit the proofs, since they are relatively simple and essentially the 
same claims, with minor differences in the coordinate system, were proved in 
[22]. (See also [3] where similar facts are proved.) 
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Lemma 8. (i) The space Qy is compact in the uniform norm. 

(a) The functional J on Q — y^Q<y<iQy is a quadratic functional which can 

he written as 

Jig) = {9,9), 

where (•, •) is defined by 

{g^h)^-l I g' {s)h' (t) log \s -t\dsdt. 
Jo Jo 

The bilinear form (■, ■) is defined for any two Lipschitz functions g,h, is 
continuous on Q with the uniform norm, and is positive semidefinite in the 
sense that {g,g) > for any Lipschitz function g. The restriction of (■, ■) to 
Qy is positive-definite. 

(iii) J is strictly convex on Qy. Therefore, a minimizer g* exists and is 
unique. 

The lemma already solves the problem in the case y = 1/2, where clearly 
gy^ = is the minimizer for J among all Lipschitz functions, and in partic- 
ular on Qi/2- It is also easy to see that a function g is the minimizer for J 
on Qy if and only if —g is the minimizer on ^i-^^. So we may assume for the 
rest of the discussion that y < 1/2. 

With these preparations, we can start the analysis. We need to minimize 
J{g) under the constraints g e Qy, which we rewrite as 

(i) ^(0) = 0, 

(ii) g is differentiable almost everywhere and g' satisfies 

-l<9'<h (17) 

(iii) lo 9'{x) dx^2y-l. 

To address the third constraint, we consider J as being defined on the larger 
space Q and form the Lagrangian 

^9, A) = ^(5) - A / g\x) dx, 
Jo 

where A is a Lagrange multiplier. Minimizing J under this constraint trans- 
lates to the equation 

W{s) := = - g'{t) log \s-t\dt-X = 0. (18) 
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This should hold whenever g'{s) is defined and is in (—1, 1). Because of 
the constraint (17), the condition will be different when g' = —1 or g' = 1. 
The correct condition (the so-called "complementary slackness" condition) 
is given by the following lemma. 

Lemma 9. If g & Gy and for some real number A the function W{s) defined 
in (18) satisfies 

r=o z/^'(s)G(-i,i), 

W{s) ts 1>Q ifg\s) = -l, (19) 
[<0 tfg\s) = l, 
then g = gy is the minimizer for J in Qy. 

Proof. We copy the proof almost verbatim from [22i Lemma 7]. \i h E Qy, 
then in particular h is 1-Lipschitz, so 

{h'{s)-g'{s))W{s)>Q 

for all s for which this is defined. So 



h'{s)W{s)ds > / g'{s)W{s)ds, 
Jo Jo 

or in other words 

2{g, h) - \{2y - 1) > 2{g, g) - \{2y - 1), 

which shows that ((7, h) > {g,g). Therefore we get, using Lemma |8](ii), that 

{h, h) = {g, g) + 2{g, h - g) + {h - g,h - g) > {g,g), 

as claimed. □ 

Having established a sufficient condition (comprised of the three separate 
conditions in (19)) for a function to be a minimizer, we first try to satisfy the 
condition (18), and save the other conditions for later. Based on intuition 
that comes from the problem's connection to the combinatorial model, we 
make the assumption that the minimizer g is piecewise smooth and satisfies 



g'{s) e (-1,1) ifsG 
g'{s) = -1 ifs^ 



1-p l + P ' 

2 ' 2 
1-/3 1 + 



(20) 
(21) 
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where 

[5 = 2^y{l-y). 

Note that g'{s) = —1 translates (via (16)) to f'{s) = in the original space 



J-'y of y-admissible functions, which corresponds to having no ascents (or 
very few ascents) in the vicinity of the scaled position (s, y) in the height 
matrix of the ASM. Our knowledge of the endpoints of the interval in which 
g'{s) > —1 is related to our foreknowledge of the arctic circle theorem, and 
one might raise the criticism that this constitutes a "guess". However, the 
analysis in [22| shows that it would be possible to complete the solution even 
without knowing this function in advance; here, we guess its value (which 
actually can be easily guessed based on empirical evidence) so as to simphfy 
the analysis slightly. 



Substituting this new knowledge about g into (18) gives the equation 



2 



i -,a 

2 



g'{t) log \s — t\ dt = X — slogs — (1 — s) log(l — s) 
+ s log — s] - (5, s G 



Differentiating with respect to s then gives 

dt = -logs + log(l - s) 



^ s — t 



2 



+ ^og[s-^—^] -\og(^-s]. (22) 



So, just like in the analysis in [22], we have reached the problem of inverting 
a Hilbert transform on a finite interval (the so-called airfoil equation). 
Moreover, the function whose inverse Hilbert transform we want to compute 
is very similar to the one that appeared in [22] - in fact, up to scaling factors 
only the signs of some of the terms are permuted, and in [22] there is an 
extra term equal to the Lagrange multiplier A. 

Now recall that in fact the general form of the solution of equations of 
this type is known. The following theorem appears in [TJ Sec. 3.2, p. 74] 
(see also [23l Sec. 9.5.2]): 
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Theorem 10. The general solution of the airfoil equation 

1 h{u) 



TT /_i U — V 



du = p{v), 



If I < 1 



with the integral understood in the principal value sense, and h satisfying a 
Holder condition, is given by 



h{v) 



1 



^ \/l — v?p{u) , c 
' du + 



TT Vl - V"^ 7-1 V - U 

for some c. 
Now set 

This function should satisfy 



(23) 



^ h{u) 



u — V 



du = log 



1-Pu 



- log 



1 + Pu 



log(l + u) - log(l - m), 



so, applying Theorem TOj we get the equation 
1 1 



h{v) 



c 

+ 



log 



1 + u 
l-u 



+ log 



l-/3u 
1 + pu 



du 



Where c is an arbitrary constant. This can be written as 

1 ..... c 



/.(f) 



:(J(f,l//?)+X(-f,l//?)) 



(24) 



where X is defined by 

^(e,7) 



and is evaluated in [221 Lemma 81 as 



log 



1+V 
7 + V 



drj 



X(^, 7) = TT 1 — 7 + a/7^ — 1 — ^arccosh(7) 

— 2a/1 — .^^ arctan 



' (7-l)(l-0 
(7 + l)(l + 
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Therefore we get that 



Kv) = r. o c + 




03-1 + 1)(1 + t.) V03-' + l)(l-") 



Since c is an arbitrary constant, we see that the only sensible choice that 
will allow to be a bounded function on the interval (—1, 1) is that of 
c = - 1 + ^1 -/52)//5. So we have 




hiv) = --{ arctan , / ^^n^yy^^ + arctan J ^^n^yy^^ 



At this point, it is worth pointing out that in (24), if we had the difference 
of the two X integrals instead of their sum, we would get at the end (up to 
some trivial scaling factors that are due to the use of different coordinate 
systems) exactly the function from the paper [22] that solves the variational 
problem for random square Young tableaux! (Compare with eq. (36) in [22] 
and subsequent formulas). Thus, while the variational problems arising from 
these two combinatorial models are not exactly isomorphic (which would be 
perhaps less surprising), they are in some sense nearly equivalent. It would 
be interesting to understand if this phenomenon has a conceptual explanation 
of some sort, but we do not see one at present. 

Simplifying the expression for h using the sum-of-arctangents identity 

X -\-Y 

arctan X + arctan Y = arctan — — 

1 — Xi 

gives 

h{v) = - 
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Going back to the original function g related to h via (23), we get that 



g'{s) = M(23-l)//3) = -^arctan./ ~\ ^\ i 



arctan 

TT 



1/2 



arctan 



TT 



7i-(y-l/2)^-(3-l/2)2^ 

V^i-(y-i/2)^-(^-i/2)^ \ _ 



for s G (^-y^, ^^)- From this, we can now get g by integration. First, from 
(21 ) we obtain that 



if < s < 



1-/3 



Next, in the interval (^-y^, ^y^) we can integrate g' using the identity 



/ arctan \J a — du = t arctan 
Jo 

+ + a arctan 



a/I + a\/a — 
and obtain without much difficulty that 

■1-/9 



— arctan 



{t^ < a). 



l-<3 
2 



y - s + 



TT 



+ - arctan 
2 



[s — 1/2) arctan 



i - (s - 1/2)2 -(y- 1/2)2 



1/2 



-(1/2 — y) arctan 



2(s-l/2)(l/2-y) 
V'i - - 1/2)^ - (y - 1/2)^ 
s - 1/2 



Vi-(^-V2)^-(i/-i/2)^ 
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for.G(Y,^)- 

Finally, from this last equation it is easy to check that 



9 



l + P 



lim g{s) = 2y- 



4^ 



so, for s > 2 5 



- again because of (21 ) we get that g{s) = 2y — s. In particular, 
g satisfies the conditions g{0) = 0, g{l) = 2y — 1, and it is also 1-Lipschitz, 
so g e Qy. 

To summarize, we have recovered as a candidate minimizer exactly the 

We also verified that it is in Qy. Furthermore, by 



function from Theorem 7 



the derivation and the use of Theorem 10, we know that it satisfies (22), or in 

f, We wanted to show that W{s] 



other words that W'(s) = on 







on this interval. But looking at the definition of W{s) in ( 18 ), we see that we 



are still free to choose the Lagrange multiplier A, which starting from eq. (22) 
has disappeared from the analysis! So, taking ^ = — Jq g'if) log \ t — \/2\dt 
ensures that (18) holds on {j^t ^^), which is one of the sufficient conditions 
in Lemma ini 

All that remains to finish the proof that g = g* is the minimizer is to 



verify the second and third conditions in (19), which we have not considered 
until now. The third condition is irrelevant, since g' is never equal to 1, 
so we need to prove that W{s), which we will now re-denote by W{s,y) to 
emphasize its dependence on y, is nonnegative when s ^ [^^) ^~2^]- Since 
g' is an even function, it follows that W{-, y) is also even, so it is enough to 
check this when s > 

Once again, our argument follows closely in the footsteps of the analogous 

and let y = ^^^2^^^ , so that 
= 0. To finish the proof, it is 



part of the proof in [22] ■ Fix 1/2 < s < 1, 

(3{y) = s. We know from (18) that W{s,y) 
enough to show that 



dW{s,y) 
dy 



< iorO<y<y. 



Denote G{x,y) = gy{x). Then 
dW{s,y) f'd'Git,y) 



dy 



dtdy 



d^G(t,y) , , 
log\s-t\dt+ I log \t 



l/2\dt. 
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A computation shows that if t G ^^-^f^, j then 



StSy dy-y--' 7, _ _ 1/2)2 _ (2,-1/2)^' 

and otherwise d'^G(t,y)/dtdy is clearly 0, so that 

dW{s,y) _2 r^^+^^/'^ log|t- 1/2| -log(s-t) 

^Al-/3)/2 - - 1/2)2 -(y- 1/2)2 

Now use the two standard integral evaluations 



log \x\ 



Mog(a-x) fa + Va^- 1\ , 

ax = vriog I , (a > Ij, 



(see pLS, eq. 4.241-7, p. 533], and [131 eq. 4.292-3, p. 553]) to conclude that 

dW{s,y) - 1/2 +^(,- 1/2)2 -(/?/2)2 

- -2 log 



5|/ " V /9/2 

Since we assumed that y < y, or in other words that s > ^'-^y^, it follows 
that 

dW(s,y) , /s-l/2\ 

^ '^^ < -2 log — — ^ < 0, 



dy - "V /3/2 
as claimed. 

We have proved Theorem 7[ (hence also Theorem [?]) , except the claim 
about the value of the integral functional J at the minimizer g*. This value 
could be computed in a relatively straightforward way, as was done for the 
analogous claim in [22j. We omit this computation, since, as was pointed 
out in [22], this can also be proved indirectly by using the large deviation 
principle to conclude that the infimum of the large deviations rate functional 
/(/) -|- 9{y) over the space J-'y must be equal to 0. Therefore the proof of 
Theorem 7' is complete. □ 
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5 The limit shape of P^Jj^j^- random ASM's 

We now apply the results from the previous sections to prove a limit shape re- 
sult for the height matrix of random ASM's chosen according to the measure 

pn 

Theorem 11. Let F{x,y) = fy{x), where for each < y < 1, f* is the 
function defined in (15). For each n let Mn be a F^^^-random ASM of order 
n, and let Hn = H{Mn) = {h'^j)^j^Q be its associated height matrix. Then as 
n —>■ oo we have the convergence in probability 



max 

0<i,j<n 



-^-F{z/n,j/n) 
n 



0. 



Proof. Fix e > 0. We want to show that 



A? 



max 

0<ij<n 



J2I 

n 



- F{i/n,j/n) 



> e 



satisfies P^omI^") ^ as n — > 00. We start by showing a weaker statement, 
namely that ii y E (0, 1) is given, then F'^ouiB'") — > as n — *• cxo, where 



B. 



max 

0<i<n 



n 



Fiy,j/n) 



> e/2 



(and [xj denotes as usual the integer part of a real number x). To prove 
this, note that 



where the union is over all (n, A;)-admissible sequences u (with k = \ny\) 
such that 



ii/u-/;i 



max 

0<x< 



:J/u(x)-/;(x)|>6/2 



(here. 



||oo denotes the supremum norm on continuous functions on [0, 1]). 
The number of such sequences is bounded by the total number of (n, k)- 
admissible sequences, so is at most (^) < 4" (since an (n, /c)-admissible se- 
quence is determined by the positions of its k ascents), and for each such u, 
by Theorem |6] we have 



PSoM {Me An-. H{M)\ny\ = u) < C exp(-c(e, y)^^ 
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where C is a universal constant, and 

c{e,y) = inf {/(/) + e{y) : f e J'y, \\f - /;|U > e/2} . 

The set oi f E J^y such that ||/ — /*||oo>e/2isa closed set in the compact 
space Qy with the uniform norm, therefore by Theorem [7] and Lemma |8} 
c(e, y) is a positive number. Combining these last observations, we see that 
indeed Pdom(-^I!?/) ^ as — > cxo. 

Next, we claim that the event is contained in the union of a finite 
number (that depends on e but not on n) of events B^^^ , so if PdomI-^I^?/) ~^ 
for all y then also ¥'^q^{A^) 0. This follows because of the Lipschitz 
property of the height matrix and of the limit shape function F, which means 
that proximity to the limit at a sufficiently dense set of values of y implies 
proximity to the limit everywhere. The details are simple, so we leave to 
the reader to check that taking yj = [je/8j for j = 1, 2, ... , [8/eJ is in fact 
sufficient to guarantee that 

^ U 

i=i 

as required. □ 

In the next section we will use a connection between uniformly random 
domino tilings of the Aztec diamond and P^g^^-random ASM's to prove a 
limit shape theorem for the height function of the random domino tiling. It 
will be helpful to consider for this purpose a variant of the height matrix of an 
ASM M, which we call the symmetrized height matrix (it is sometimes 
referred to as the skewed summation of M). If M G An-, we define this 
as the matrix Hsym{M) = {h*^^)'^^^^ with entries given by 

hlj='+3 - ^H{M)i,,, {Me An, 0<t,j<n), 

where H{M)ij is the {i,j)-th entry of the (ordinary) height matrix of M. 
See Figure [Tol for an example. The following theorem is an equivalent version 
of Theorem [TT] formulated for these matrices. 
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/0 1 23456\ 

1 2 3 2 3 4 5 

2 3 2 3 4 3 4 

3 2 3 4 3 2 3 

4 3 2 3 2 3 2 

5 4 3 2 1 2 1 
\6543210/ 

Figure 10: The symmetrized height matrix of the ASM from Figure [7j 



Theorem 11'. Let G{x,y) = x + y — 2F{x,y), where F is defined in The- 



orem 

m = 



11 



For each n let Mn be a F^Q^-random ASM of order n, and let 
HsYuiMn) = (/i*j")"j=o be its associated symmetrized height matrix. 
Then as n ^ oo we have the convergence in probability 



max 

0<i,j<n 



h* 



n 



G{i/n,j/n) 



0. 



We remark that it would have been possible to work with symmetrized 
height matrices right from the beginning. In that case the large deviation 
analysis would have lead directly to Variational Problem[2]without going first 
through Variational Problem [l] (Note that the limiting symmetrized height 



function G{x,y) can also be written as G{x,y) 
solution to Variational Problem U\ 



y ~ 9l{^)^ where g* is the 



6 Back to domino tilings 

We now recall some basic facts from [B] about domino tilings of the Aztec 
diamond AD„, their height functions, and their connection to alternating sign 
matrices and their height matrices. This will enable us to use our previous 
results to reprove the Cohn-Elkies-Propp limit shape result for the height 
function of a uniformly random domino tiling of AD„ as n — * oo. 
Let Q = Q{ADn) be the directed graph whose vertex set is 

r(AD„) = {(2,j)eZ2:|z| + |j|<n + l}, 
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and where the adjacency relations are 

(hJi) («2,j2) or 

ii = i2 and ji - = 

We call Q{KDri) the Aztec diamond graph. Note that its adjacency 
structure is the standard nearest-neighbor graph structure induced from 
Z^, where in addition edges are directed according to a checkerboard par- 
ity rule, namely, that if a checkerboard coloring is imposed on the squares 
[n, n -|- 1] X [m,m + 1] in the lattice dual to Z^, then the nearest-neighbor 
edges u — > f are all directed such that a traveller crossing the directed edge 
will see a black square on her left; see Figure 11 ^a). 

Define a height function to be any function rj on V(AD„) such that for 
any edge u ^ v m. ^(AD„) we have 

Tjiu) — Tjiy) = 1 or — 3, 

and such that ri{u) —rjiv) = 1 whenever u ^ v is one of the boundary edges. 
A height function t] on y(AD„) is called normalized if ri{—n, 0) = 0. 

It is known that any domino tiling T of AD„ determines a unique normal- 
ized height function 7]t by the requirement that for any directed edge u ^ v 
we have 



riTKu) - r]T[v) 



-3 the segment (u,v) crosses a domino tile in T, 
1 otherwise. 



Conversely, any normalized height function rj is of the form rjT for some 
domino tiling. See Figure [Tl|b). 

Another important fact concerns the beautiful connection, discovered by 
Elkies, Kuperberg, Larsen and Propp [6] , between height functions of domino 
tilings of AD„ and height matrices of ASM's: each normalized height func- 
tion T] on V^(AD„) is essentially comprised of the superposition of two (sym- 
metrized) height matrices Hsym{A) , Hsym{B) where A is an ASM of order n 
and B is an ASM of order n + 1. More precisely, Hsym{A) and Hsym{B) can 
be recovered from r] by 

HsUB)., = + (26) 



38 



8 



6 — 7 — 6 
4 — 5 — 4 — 5 — 4 
2 — 3 — 2 — 3 — 2 — 3 — 2 
— i — 4 — 5 — 4 — i — 4 — i — 
2 — 3 — 6 — 3 — 2 — 3 — 2 
4 — 5 — 4 — 5 — 4 
6 — 7 — 6 
8 

(b) 

Figure 11: (a) The Aztec diamond graph of order 3; (b) The normahzed 
height function corresponding to the tihng from Figure [T} 

(note the shght difference from the formulas in [6] due to a difference in the 
center of the coordinate system used) . This correspondence defines a one-to- 
one mapping from the set of domino tihngs of AD„ to the set of pairs {A, B) 
where A G An and B G An+i- The pairs {A, B) which are obtained via this 
mapping are exactly the so-called compatible pairs defined by Robbins and 
Rumsey [25] : A and B are called compatible if the (non-symmetrized) height 
matrices H{A),H{B) satisfy the conditions 

HiB),, < H{A\,, 

if(5),+i,,+i - 1 < H{A\,, 

H{A),j < H{B),+i,,, 

H{A),j < H{B\,j+i. 

It was also shown in [25] that for a given ASM A G An-, the number of 
B G A^n+i that are compatible with A is equal to 2^+^"^^ Combined with the 
formula for the number of domino tilings of AD^, this implies that if T is a 
uniformly random domino tiling of AD„, and {A, B) is the associated pair 
of compatible ASM's, then the random ASM A is distributed according to 
the domino measure Pdom (of course, this provides the explanation for our 
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choice of name for this measure). 



We now combine Theorem |11[| with the above discussion to easily obtain 
the following result, originally proved in [2]. 

Theorem 12. For each n > 1, let Tn be a uniformly random domino tiling 
of ADn, and let rjn = tjt^ be its associated height function. Then as n oo 
we have the convergence in probability 



max 

ii,j)eViAD„) 



-Vn{i,j) - R{i/n,j/n) 



n 



where 



R{u,v) = 2G 



u — V 



1 u 



{\u\ + \v\<l), 



and G is defined in Theorem 11 



Proof. For pairs (i, j) G V(AD„) for which i + j + n is odd, the proximity 
of n~^rin{i,j) to R{i/n,j/n) follows from (25). For other pairs {i,j), apply 
the previous observation to an pair {i',j') adjacent to {i,j) and use the facts 
that \rin{i,j) — ^ 3 and that i? is a continuous function. □ 



7 Concluding remarks 

7.1 Relation to the arctic circle theorem 

Theorem 12 implies a weak form of the arctic circle theorem (Theorem [T|: 
First, since inside the arctic circle the limit shape function R{u,v) is not a 
linear function, it follows that the frozen region cannot extend in the limit 
into the arctic circle, which is "half" of the theorem. In the other direction, 
we get only a weaker statement that outside the arctic circle we can have 
in the limit at most o(n^) "non-frozen" dominoes, since that is what the 
linearity of the limiting height function in that region implies. 

It is interesting to contrast this with the square Young tableaux problem. 
There, too, the large deviation approach gave only a bound in one direction 
on the behavior of the square Young tableau along the boundary of the 
square. However, Pittel and Romik managed to prove the other direction 
using an additional combinatorial argument (inspired by a method of Vershik 
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and Kerov [21]). It would be interesting to see whether one can emulate this 
approach in the present case to get a new proof of the arctic circle theorem. 
A similar question applies to the problem of random boxed plane partitions 
studied by Cohn, Larsen and Propp |3j , where again the limit shape theorem 
for the height function does not imply an arctic circle result in its strong 
form. 

7.2 Other arctic circles 

In this paper we have shown that two so-called arctic circle phenomena, 
namely those appearing in the contexts of random domino tilings of the 
Aztec diamond and of random square Young tableaux, are closely related, 
in the sense that the limit shape results underlying them can be given a 
more or less unified treatment using the techniques of large deviation theory 
and the calculus of variations, and that the derivations in both cases result in 
nearly identical computations and formulas. Note that these are not the only 
combinatorial models in which arctic circles appear. Other examples known 
to the author include the shape of a uniformly random boxed plane partition 
derived by Cohn, Larsen and Propp [3] and the arctic circle theorem for 
random groves, due to Petersen and Speyer [21]. One might therefore wish to 
extend the insights of the present paper to these other models. The treatment 
of boxed plane partitions in |3] is already based on a large deviations analysis, 
and in fact the variational problem studied there seems to be quite closely 
related to the variational problems studied here and in |22]. Therefore, it 
should be relatively straightforward to use the techniques presented here 
to give a new derivation of the solution to the variational problem from 
[3] (which in particular would provide a fully satisfactory answer to Open 
Question 6.3 from that paper). 

The analysis of random groves, on the other hand, is based on generating 
function techniques, and it is not clear how to apply the ideas presented here 
to that setting. 

7.3 Uniformly random ASM's 

One reason why the methods and ideas presented in this paper may be con- 
sidered worthy of attention is somewhat speculative in nature. It pertains to 
the potential future applicability of these methods and ideas to a well-known 
open problem on alternating sign matrices: that is, the problem of finding 
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Figure 12: The Colomo-Pronko conjectured limit shape for uniformly random 
alternating sign matrices. 



the limiting shape of a uniformly random ASM of high order. Here, "limit 
shape" is usually taken to refer to the shape of the region in which the non- 
zero entries cluster (the "temperate region"), although one could also ask (as 
we have done here in the case of P^Q^j-random ASM's) about the limiting 
shape of the height matrix, which also contains useful information about the 
behavior of the ASM inside the temperate region. 

Important progress on this question was made recently by Colomo and 
Pronko [5], who conjectured the explicit formula 

+ + \xy\ = \x\ + \y\ 

for the limit shape of the boundary of the temperate region in a uniformly 



random ASM (Fig. 12), and provided a heuristic derivation of this conjec- 
tured formula based on certain natural, but still conjectural, analytic as- 
sumptions. 

In view of this state of affairs, it is worth noting that the ideas presented 
in this paper seem to be rather suitable for attacking this challenging open 
problem. There is only one main "missing piece" (albeit possibly a very 
substantial one) in our understanding. The idea is to replace Theorem |3} 
which is the combinatorial observation which lies at the heart of the large 
deviations analysis, with an analogous statement that holds for the uniform 
measure on the set An of ASM's of order n. This statement is given in the 
following theorem, whose proof follows similar lines to the proof of Theorem [3] 
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and is omitted. 



Theorem 13. Let Punif denote the uniform measure on the set of ASM's of 
order n. For a positive integer k and integers Xi < X2 < ■ ■ ■ < Xk, denote by 
ak{xi, . . . ,Xk) the number of monotone triangles of order k with bottom row 
(xi, . . . ,Xk)- Then, in the notation of Theorem^ we have 



Unif 



MeAn-. (Xfc(l), . . . , Xk{k)) = (xi, ...,Xk) 

1 



\A I 



. . . , Xfc) On-kiVl, • • • , Vn-k) 



Unfortunately, while a formula for \An\ is known (see pj), the function 
ak seems much more difficult to understand (and in particular, to derive 
asymptotics for) than the Vandermonde function A, and this is the piece 
that is missing when one tries to duplicate our analysis to the setting of 
uniformly random ASM's. Nevertheless, the function has recently been 
the subject of several very fruitful studies. Fischer |9j derived the following 
beautiful "operator formula" for ak- 



ak{xi, ...,Xk) 



l<i<j<k 



A{xi,...,Xk) , . 

A{l,...,k) ■ ^'^^ 



Here, A is the Vandermonde function as before, and Id, Ej and Di are opera- 
tors acting on the ring of polynomials C[xi, . . . , Xk]'- Id is the identity opera- 
tor, Ej is the shift operator in the variable Xj (that substitutes Xj + 1 for each 
occurrence Xj in a polynomial), and Di = Ei — Id is the (right-) differencing 
operator in the variable Xi. 

Fischer then showed in several subsequent papers that it is possible to use 



(27) to get highly non-trivial information on the enumeration of alternating 
sign matrices: In [TU] she obtained a new proof of the celebrated Refined 
Alternating Sign Matrix Theorem (see [1] for the statement and fascinating 
history of this result); in [12] she and the author obtained additional results 
concerning a "doubly-refined" enumeration of ASM's; and in [TT] she proved 
a conjecture from [12], extending these results further to a "multiply-refined" 
enumeration. Thus, it seems quite conceivable that additional study of 
may eventually lead to a deeper understanding of this function, that, in 
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combination with Theorem 13 and the techniques of this paper, could provide 
a basis for a successful attack on the limit shape problem for uniformly 
random ASM's. 
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